Inverse scattering problem for quantum graph vertices 
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We demonstrate how the inverse scattering problem of a quantum star graph can be solved by 
means of diagonalization of Hermitian unitary matrix when the vertex coupling is of the scale 
invariant (or Fiilop-Tsutsui) form. This enables the construction of quantum graphs with desired 
properties in a tailor-made fashion. The procedure is illustrated on the example of quantum vertices 
with equal transmission probabilities. 
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The interest in the inverse scattering problem for quan- 
tum graphs [1-3 is two-fold. These graphs are a prime 
example of solvable systems possessing nontrivial physi- 
cal properties [4 . At the same time, it is also important 
for its relevance as the design principle of nanowire-based 
single electron devices. 

In this letter, we consider the inverse scattering prob- 
lem on a star graph with scale invariant vertex coupling 
[5] , which is an important subset among all the couplings 
preserving the probability current [6 ; we note that the 
corresponding scattering matrix is energy independent. 
A star graph with a number of half-line edges connected 
at a single point is the elementary building block of a 
generic graph. In general, the inverse problem for star 
graphs is easy to solve [7 but the solution does not give 
much insight into the physical meaning of the coupling. 

Here we present an alternative approach for star graphs 
with scale invariant coupling. We shall give the solution 
to the corresponding inverse scattering problem in the 
form of eigenvalue problem of a Hermitian unitary ma- 
trix. In particular, we consider quantum vertices with 
equal transmission probabilities, a subclass of the scale 
invariant case, and we show that the problem boils down 
to a search for a special class of Hermitian unitary matri- 
ces which can be regarded as a generalization of complex 
Hadamard matrices. Intriguing designs emerge for the 
realization of quantum devices with such properties. 

Note that any quantum-graph vertex coupling behaves 
effectively as a scale-invariant one in both the high and 
low energy limits [8 , which gives a hint that our analysis 
might have an extension to the general case. 

Consider thus a star graph vertex of degree n, with n 
half-line edges sticking out of a point-like node. The scale 
invariant subfamily of the general coupling conditions is 
characterized by a complex matrix T of size (n — m) x m, 
where m can take an integer value m = l,2,...,n — 1, 
being given by 
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where I^^^ signifies the identity matrix of size / x and 



the boundary- value vectors ^ and ^' are defined by 
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with il)i{xi) and il^[{xi) being the wave function and its 
derivative on i-th edge [9]. The coordinates Xi on the i- 
th half-line are labeled outwardly from the vertex, which 
corresponds to = for all i. To cast the coupling into 
the form ([T]) one has, in general, to renumber suitably 
the edges. The quantum particle coming in from the j- 
th edge and scattered off the vertex is described by the 
scattering wave function on the i-ih line, ip\'^\x)^ which 
is of the form 

^/;p')(x) = 5,,,e-"="+5,je"=^ (3) 
Let us express S in terms of T. From ([3| we have 

^(0) = / + 5 and ^'(0) = k{-I + 5); 
we substitute into ([T]), which leads to the equation 
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It is easy to observe from Q that 

with the matrices X^, defined by 

_ T \ \ 
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We see that ([5| can be viewed as a diagonalization 
formula of Hermitian unitary matrix S with a diagonal- 
izing matrix of specific block diagonal form, X^, which 
gives a prescription to obtain T that defines the boundary 
condition from the scattering matrix S. In other words, 
solution of the our inverse scattering problem is given in 
terms of a diagonalization. 

In practice, the procedure of obtaining T by a diago- 
nalization of S can be cumbersome for large n, and there 
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is an alternative simpler way. A calculation shows that 
([5| can be rewritten in the form 
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Let us divide S into four submatrices 5i2, ^21 and 
of size m X m X {n — m)^ {n — m) x m and (n — 
m) X {n — m)^ respectively, as follows 



S 



S21 S22) 



(8) 



The block matrices Sij express in terms of T as 5ii = 

-/^^) + 2 (/^^) +TTt)"\ S12 = 2 (/^^) +TTt)"^T, 

and 522 = /^^-^) - 2 (/(^-^) +TtT)"\ From here, 
one gets easily 

T = + 5n) "'^12= Sl^ (/("-™) - 522) . (9) 

Hence, the algorithm to obtain the matrix T character- 
izing the vertex is the following: 

1. Take the scattering matrix S and set m = rank(tS + 

2. Decompose S according to ([8|. If necessary, change 
the numbering of the incoming edges so that /^^^ + 
Sii is regular (it is always possible). 

3. Calculate T using 

We remark that the matrix T obtained by the algorithm 
above naturally depends on the numbering of the edges 
we choose. 

In the rest of the paper we demonstrate how the matrix 
T is used for understanding the meaning of the coupling 
and for the construction of the vertex with prescribed 
scattering properties. 

Asking about the meaning of the coupling, recall a star- 
shaped network with a potential at the node can tend 
in the zero-diameter limit to the star graph with S cou- 
pling [To]. More complicated couplings can be obtained 
by using 6 vertices as building blocks, applying localized 
magnetic fields to achieve phase change if necessary; the 
aim is to devise a design principle to construct an ar- 
bitrary coupling condition. Here the matrix T derived 
above plays an important role. In the scale-invariant 
case ([1]) with the elements of T = {t^j}, i = l,...,m 
and j = m + 1, ..,n, given, the scheme works as follows 

(i) Take the endpoints of the n edges, numbered by 
j = 1, 2, n, and connect them in pairs (z, j) by internal 
edges of length d/r^j, except when r^j = in which case 
the pair remains unconnected. Apply a vector potential 
Aij on the segment (i, j) to produce extra phase shift Xij 
between the endpoints when its value is nonzero. Place 
the S potential of strength Vi at each endpoint i. 



(ii) The length ratio r^j used above and the phase shift 
Xij are determined from the non-diagonal elements of the 
matrix Q defined by 
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using the relation r^^e^^^^ = Qij {i ^ j). This means that 



we have Tije^^'^ 



tij for i < m^j > m, and r^je^^ 



Xl^>m ^'ii'^'ji — ^5 i,j>m we have r^j = and 

naturally also Xij = 0- 

(iii) The S coupling strength Vi is given by the diagonal 
elements of the matrix V defined by 



1 



(11) 



where R is the matrix whose elements equal the absolute 
values of the matrix elements of Q, i.e. R = {vij} — 
{IQijI}; the nxn matrix J^^^ has all the elements equal 
to one. This means that we have '^^ = ^(l — X^Km 

for i > m, and = ^ ( Y.i>mVl - - Eii^i)<rJii) for 
i < m. The described way to fine tuning of the lengths 
and 5 coupling strengths is devised to counter the generic 
opaqueness brought in with every addition of a vertex or 
a connecting edge into the graph. 

The wave function 0(x)= (j)ij{x) on any internal edge 
with indices (i^j) has to satisfy the relation 
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with F{x) = xcotx and G{x) = xcosecx. Combining 
(12) with the condition at the i-th endpoint where we 



have the ^-potential of strength v^. 



(13) 



we obtain the relations between the boundary values ipi 
?/^i(0) and ip[ = '0^(0) in the form 




'ruGii^lJu (14) 



where the obvious notations Fij 



cot ^ and Ga = 

ru Til V 

— cosec — have been adopted. Note that the equation 
( p3| ) is exact and does not involve any approximation. 
For small values of the length parameter d we have Fij = 
1 + 0{d?) and Gij = 1 + 0{(P); then we can show by a 
straightforward computation in the manner of [9 that the 
shrinking limit d ^ gives the desired coupling condition 
for the scale-invariant vertex ([T]). 

With the procedure described above, it is possible to 
construct a star-graph from any given scattering matrix 
of the considered class. Our previous result detailed in 
pT] which provides a reconstruction of the "free-like" 
scattering is one such example, and it could be achieved 
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by 



FIG. 1. Finite approximation to the reflect ionless scale in- 
variant vertices corresponding to ( |16| (left) and (18) (right) 
constructed according to (lO)-(ll). Dotted line indicates the 



presence of a non-zero phase shift Xij- 



more easily by the current method. Here, we will illus- 
trate the application of the procedure on the following 
exemplary problem. Let us look at this question: Can 
one construct a quantum vertex for which the particle 
incoming from any line is transmitted to all other lines 
with equal probability? At first, we should ask about the 
existence condition for the scattering matrix of the form 



S 



1 



+ n - 1 



, (15) 



-d J 



with a non-negative real parameter d. 

These matrices have been recently examined in fi3\ . It 
has been proved that the parameter d is always bounded 
from above by | — 1 (except for n < 2), and moreover, 
that for most values of d G [0, ^ — 1], the existence of the 
corresponding matrix ([T5| is impossible if the order n is 
odd. By contrast, if n is even, then one can construct an 
S for infinitely many values of d, in particular for any d 
from the interval [| — 3, | — 1], or, under certain extra 
condition, for all d G [0, 1]; for further details and explicit 
matrix constructions we refer to p!3|. 

With a matrix of the type ([l5| in hand, we can proceed 
to the construction of the finite approximation. We will 
demonstrate it on two examples. We first look at the case 
of a reflectionless scattering with uniform transmission to 
all the other edges. In other words, we ask whether there 



is a graph with the scattering matrix (15) with d = 0. 



Vertices yielding such scattering matrices are known to 
be useful in investigation of semiclassical properties of 
quantum-graph spectra [12]. Let us adopt a result from 
[13] which says that such an S can exist only for even 
n. We limit ourselves for simplicity to real <S; we remark 



that a real matrix of the type (15) with zero diagonal 



is called symmetric conference matrix^ and is known to 
exist for n = 2, 6, 10, 14, 18, 26, 30, 38, ... . 

Here we inspect the example of n = 6 when S is given 



1 / /(3) _ J(3) _2/(3) + J(3) 



(16) 



Applying (|9]), the corresponding T is easily calculated, 

T = -7/(^) + (1 + 7)/^) (17) 

where 7 = (a/5 — l)/2 is the golden mean. Our finite 
approximation is specified by following parameters; ri2 = 
^23 = = 4 + 37,ri4 r25 rse = l,ri5 
r26 ^24 rsi = r32 = 1 + 7, r^B = r^e r^Q = 
0,6^^12 = gix23 = gixis = gixii = I for all others, and 

vi=V2=v^ = -6^, ^4 = ^5 = ^6 = -2^. The 
finite graph approximation is schematically illustrated in 
the left side of Figure 1. 

Second example is the equal-scattering graph, in which 
the scattering is uniform to all the edges including the 
one of the incoming particle. Such a matrix, called 
symmetric Hadamard matrix^ is known to exist for n = 
2, 4, 8, 12, 16, 20, 24, ... . An example of such 5 forn = 8 
is given by 

1 /2/W-j(4) _2/(4) + j(4)\ 

- 71 + _2/(4) + j(4) j • y^^) 

The matrix T specifying the vertex coupling is found to 
be equal to 



T = 



/(4) + ^ J(4) 



Cr- 1 



a + 1 cr + 1 



(19) 



where a = y/2 — 1 is the silver mean. Our finite 
approximation is specified by following parameter val- 
ues; ri2 = ri3 = ri4 = r23 = ^24 = r34 = 1 + a. 



^15 
^28 



r26 

^25 



^48 
^35 
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= ^45 
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= ^46 



ri8 = r27 
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1X12 — 



r47 
= 0, e 



olX46 — 



-1, 



^3=^4 = -^, 



^56 ^57 ^58 r67 ^68 ^78 

giXis — giXl4 — giX23 — giX35 — gi028 — g^. 

gixii — I oil others, and vi = V2 

= vq = vj = = The finite graph approx- 

imation is schematically illustrated in the right side of 
Figure 1. 

Finally, let us discuss the convergence of the the de- 
scribed finite-size graph approximation. In Figure 2, we 
display scattering matrix elements of the finite graph con- 
structed to approximate the equal-scattering reflection- 



less matrix (16). They are calculated directly from (14). 



The scale of the wave length k is given by 1/d. The ap- 
proximation can be seen to be quite good below kd < 0.2. 
Numerical analysis of other examples of different graphs 
gives essentially the same conclusion, namely that the de- 
scribed construction does represent a physical realization 
of scale- invariant vertex couplings. 

Thus the problem of finding the desired property of 
scale invariance is turned into a mathematical question 
about a Hermitian unitary matrix, and the search for sys- 
tems with S having interesting specifications other than 
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FIG. 2. Scattering probabilities as functions of incoming 
momentum k (in the unit of l/d) of finite quantum graph ap- 
proximating the equal-transmitting reflectionless vertex rep- 
resented in the left side of Figure 1. 

those examined here should follow. Also, a study of the 
bound state spectra is one thing we have completely ne- 
glected in this letter; applications to non-quantum waves, 
including particularly electromagnetic and water waves 
should be another interesting subject. 

In our finite approximation of star graph with no inter- 
nal edges, we have actually studied the low-energy prop- 
erties of graphs with internal edges all of which are con- 
nected to the external ones, which we might term depth- 
one graphs. The examination of depth-two graphs and 
beyond seems to be a natural future direction. Our re- 
sult showing the full solution to the inverse scattering 
problem is, in a sense, a partial fulfillment of the hope 
that quantum graph somehow could be a solvable model 



and useful design tool at the same time. 
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